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Abstract 

We present an X-ray anaiysis on the Gaiactic supernova remnant (SNR) G298.6-0.0 with 
Suzaku. The X-ray image shows a center-fiiied structure inside the radio sheii, impiying this 
SNR is categorized as a mixed-morphoiogy (MM) SNR. The spectrum is weii reproduced by a 
singie temperature piasma modei in ionization equiiibrium, with a temperature of 0.78 (0.70- 
0.87) keV. The totai piasma mass of 30 M© indicates that the piasma has intersteiiar medium 
origin. The association with a GeV gamma-ray source 3FGL J1214.0-6236 on the sheii of the 
SNR is discussed, in comparison with other MM SNRs with GeV gamma-ray associations, it 
is found that the fiux ratio between absorption-corrected thermai X-rays and GeV gamma-rays 
decreases as the MM SNRs evoive to iarger physicai sizes. The absorption-corrected X-ray 
fiux of G298.6-0.0 and the GeV gamma-ray fiux of 3FGL J1214.0-6236 cioseiy foiiow this 
trend, impiying that 3FGL J1214.0-6236 is iikeiy to be the GeV counterpart of G298.6-0.0. 
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1 Introduction 

The origin of Galactic cosmic rays has remained an unresolved question since their discovery. One 
of the most plausible acceleration sites is at the shock front of supernova remnants (SNRs). The 
discovery of synchrotron X-rays (Koyama et al. 1995) and very high energy (VHE) gamma-rays 
(Aharonian et al. 2004) from shells of SNRs gave some of the first clear evidence of active particle 
acceleration in these shocks. The small-scale structure and time variability of synchrotron X-rays 
imply that the acceleration is very efficient (Bamba et al. 2003; Bamba et al. 2005; Uchiyama et al. 
2007). However, it remains an open question how accelerated particles escape from the shocks to 
become cosmic rays. 

Recently, Fermi detected GeV gamma-rays from several SNRs (Abdo et al. 2010a; Abdo et 
al. 2010b; Abdo et al. 2010c). Interestingly, many of the GeV-detected SNRs are many times older 
than the non-thermal X-ray-detected or VHE-gamma-ray-detected SNRs (Funk 2011). Moreover, 
their GeV spectra are soft and show a spectral break at energies around 1-10 GeV, implying that high 
energy particles have already lost their energy or have escaped from the acceleration sites. If the latter 
is true, their soft spectra show the evidence of particle escape. However, a larger sample of GeV 
gamma-ray SNRs is needed to judge the particle escape scenario. 

We searched for GeV gamma-ray sources associated with cataloged Galactic SNRs in the 
Fermi Farge Area Telescope Third Source Catalog (The Fermi-FAT Collaboration 2015). There are 
60 GeV gamma-ray sources which spatially coincide with the Galactic SNRs, in addition to the 12 
SNRs identified by extended gamma-ray emission. 3FGF J1214.0—6236 is associated with the radio 
SNR G298.6—0.0, detected at 408 MHz and 843 MHz with a flat radio spectral index of —0.3 (Shaver 
& Goss 1970; Kesteven & Caswell 1987; Whiteoak & Green 1996). Reach et al. (2006) reported a 
possible detection of infrared emission from the direction of G298.6—0.0, implying that the shock 
of this remnant may have encountered a high-density medium. These characteristics are common 
among GeV gamma-ray-emitting SNRs. However, no detection has been made in the X-ray band 
with ROSAT PS PC (Hwang & Markert 1994). X-ray information is critical to estimate the physical 
parameters of the SNR, such as age, explosion energy, and the amount of accelerated electrons. Thus, 
our aim is to identify X-rays from G298.6—0.0 and measure parameters using the low and stable 
background X-ray observations provided by Suzaku (Mitsuda et al. 2007). 
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In this paper, we report the first X-ray imaging speetroseopy of G298.6—0.0 using Suzaku. 
The observation details are summarized in §2. An analysis of the Suzaku X-ray data is presented in 
§3, the results of whieh are diseussed in §4. 

2 Observations and Data Reduction 

Suzaku observed G298.6—0.0 on 2012, Aug. 11-13 (sequenee number: 507037010) and 2013, Feb. 
18-19 (sequenee number: 507037020). Only three set of the onboard X-ray Imaging Speetrometers 
(XISO, I, 3; Koyama et al. 2007) eould be used in this paper, due to a known problem with XIS2. 
XIS1 is a baek-illuminated CCD, whereas the others are front-illuminated. The XIS instruments were 
operated in the normal-eloeking full-frame mode. Spaeed-row eharge injeetion teehnique and the 
relevant energy ealibration method were applied to mitigate the long-term degradation of the energy 
resolution (Uehiyama et al. 2009). Sinee the observing mode is the same for the two observations, we 
analyzed the eombined data using HEADAS software version 6.16 and XSPEC version 12.8.2. We 
reproeessed the data set with the ealibration database version 2012-11-06 for XIS and 2011-06-30 for 
X-ray Teleseope (XRT; Serlemitsos et al. 2007). 

In the XIS data sereening, we filtered out data aequired during passages through the South 
Atlantie Anomaly (SAA), with elevation angles with respeet to the Earth’s dark limb below 5 deg, 
or with elevation angle to the bright limb below 20 deg in order to avoid eontamination by emission 
from the bright limb. The remaining exposure is 57 ks (17 ks for the former observation and 40 ks for 
the latter observation). 

3 Results 

3.1 Images 

Eigure 1 (a) shows 843 MHz (Whiteoak & Green 1996) and the XIS I-I-3 0.5-5 keV band images 
of G298.6—0.0 region. XIS 0 image is not used due to the dead eolumn. The position of the Eermi 
souree, 3EGE JI214.0—6236, is also shown. The radio emission shows the elear shell-like strueture, 
while X-ray emission has a eenter-filled morphology inside the radio shell and does not show any 
prominent point sourees. These results show that G298.6—0.0 is either a SNR with a eentral pulsar or 
a mixed-morphology SNR (MM SNR; Rho & Petre 1998). The point souree loeated in the north-east 
of the SNR is eoineident with the ROSAT souree IRXS J121248.7—623027 (Euhrmeister & Sehmitt 
2003) and eategorized as a rotationally variable star (Kiraga 2012), whose details are out of the seope 
of this paper. 
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3.2 Spectra 


The source spectra of G298.6—0.0 were extracted from circular region with the radius of 3.9 arcmin, 
whereas the background spectra were taken from the source free region, as shown in Figure 1 (b). We 
added XISO and XIS3 spectra for the better statistics, whereas XISl spectrum was treated separately 
due to the different response from the others. 

Background-subtracted spectra are shown in Figure 2. The spectra are rather soft and have 
emission lines, implying that the emission has a thermal origin. We thus adopted an optically thin 
thermal plasma model in collisional ionization equilibrium (VAPEC) affected by interstellar absorption 
(PHABS; Balucinska-Church & McCammon 1992). For both the emission and absorption models, the 
solar abundance values of Anders & Grevesse (1989) were used. The first fit with fixed abundances at 
the solar values was rejected with the reduced of 228.6/117 and wavy residuals. We thus fit with the 

abundances of Mg, Si, S, and Fe as free parameters. This improved fit gave a reduced of 136.1/113. 
The best-fit models and parameters are shown in Figure 2 and Table 1, respectively. We also tried to 
fit with non-equilibrium ionization plasma models, VNEI for ionizing plasma and VVRNEI for recom¬ 
bining plasma. The abundances of Mg, Si, S, and Fe were treated as free parameters. For a VVRNEI 
model, we fixed the initial plasma temperature to be 3 keV. The fitting showed similar reduced 
of 136.5/112 (VNEI) and 131.5/112 (VVRNEI). The best-fit values of the ionization/recombination 
timescale are > 9.8 x 10^^ cm“^s for the VNEI model and 6.2 (5.3-14) xlO^^ cm“^s for the VVRNEI 
model, respectively. The latter gives a slightly shorter relaxation timescale than that in ionization 
equilibrium (10^^“^^ s cm“^). However, the improvement of the fit is only marginal; both models 
show similar reduced x^- Moreover, there is no clear indication of strong radiative recombination 
continuum in the residual which is the direct evidence of a recombining plasma (Yamaguchi et al. 
2009). Thus we concluded that the plasma is well described by an absorbed plasma emission model 
in ionization equilibrium. We note that the XIS1 spectrum has residuals in the low energy band, which 
may be due to the calibration uncertainty of contamination. 

4 Discussion 

4.1 Physical Parameters of G298.6-0.0 

We report the first detection of thermal X-rays from inside the radio SNR G298.6—0.0. The charac¬ 
teristics of G298.6—0.0 classify it as a new MM SNR. Here, we estimate several physical parameters 
for this SNR. 

First, we estimate the distance to this source. We measured the absorption column to the 
source to be 1.65 (1.45 — 1.99) x 10^^ cm“^. The total hydrogen density throughout the Galaxy in 


this direction is estimated by the method of Willingale et al. (2013). The Hi column density in this 
direction is 1.4 xlO^^ atoms cm“^ (Kalberla et al. 2005), whereas the column density of molecular 
hydrogen (H 2 ) is 1.4x 10^^ atoms cm“^ from the dust map (Schlegel et al. 1998). The total absorption 
column is thus 1.5x10^^ atoms cm“^, which is similar value to our target. This indicates that the 
SNR is located in the outer part of the Galaxy. According to the face-on map by Nakanishi & Sofue 
(2006), most of the interstellar hydrogens in this direction are located within 10 kpc from the Sun. 
Thus we adopt the distance of 10 kpc. This is an independent constraint from that indicated by the 
Tj — D relation in the radio band (9.5 kpc; Case & Bhattacharya 1998). The previous method only 
gives a rough estimate and more detailed studies are necessary to get a better measurement of the 
distance to SNR G298.6—0.0. 

Assuming the plasma fills a 3.9 arcmin sphere, the total volume is 1.8xl0^®i9fQ cm^, where 
Diq is the distance in the unit of 10 kpc. The average density derived from the emission measure 
(Table 1) is cm“^, where / is the filling factor. The total mass and thermal energy 

can be estimated as and 6.4 x erg, respectively. The large total mass 

indicates that the main component of the plasma originates from swept-up interstellar medium. 

4.2 Comparison with other SNRs with GeV counterpart 

Here, we discuss on the comparison with G298.6—0.0 and other MM SNRs with GeV gamma rays, 
to judge whether 3FGL J1214—6236 has similar properties to other GeV counterparts of MM SNRs. 

GeV emission is coincident with the shell of G298.6—0.0, which is same to other SNR cases 
(Abdo et al. 2010a). The possible interaction with a high-density medium from the direction of this 
SNR is reported (Reach et al. 2006), which is also similar to other GeV SNR samples. On the other 
hand, we have no information which part of this SNR interacts with high-density medium. Further 
studies in infrared or molecular cloud bands is needed to conclude the interaction. The GeV spectrum 
of 3FGL J1214.0—6236 has a spectral break around a few GeV (The Fermi-LAT Collaboration 2015), 
which is also similar to the spectra of other GeV SNRs. These facts imply that 3FGL J1214.0—6236 
is a possible GeV counterpart of G298.6—0.0. 

Table 2 shows the summary of the properties across the absorption-corrected 0.3-10 keV and 
0.1-100 GeV bands of the GeV gamma-ray emitting MM SNRs. Large distance uncertainties for 
these sources make it difficult to apply a detailed model of the X-ray and GeV gamma-ray evolution 
to the case of G298.6—0.0. Instead we compare the flux ratio between absorption-corrected X-rays 
and GeV gamma-rays, which is independent from the distance. 

We selected samples with known total absorption-corrected X-ray and GeV gamma-ray fluxes. 
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Recent studies revealed that some of our samples have X-ray emission from recombining plasma with 
detailed analysis (Ozawa et al. 2009; Uchida et al. 2012; Sato et al. 2014), but we concentrate on total 
X-ray flux from older results as it is an easy comparison for the other fainter sources. IC443 is a 
famous MM SNRs with GeV gamma-rays (Yamaguchi et al. 2009; Abdo et al. 2010c), but do not 
have a good estimate of the total X-ray flux in the literatures due to their large apparent diameters, 
and thus we did not adopt it as our samples. We also ignored emission from ejecta reported in W49B 
and G349.7-I-0.2 (Kawasaki et al. 2005; Yasumi et al. 2014), which contributes only ~10% flux of the 
total flux. GeV gamma-ray spectra for several sources have spectral breaks around 10 GeV (Abdo et 
al. 2010a; Abdo et al. 2010b; Abdo et al. 2010c, for example). To treat fainter sources equally, we 
used the total energy flux from the Fermi source catalog (The Fermi-LAT Collaboration 2015). The 
difference in the fluxes of W44 between the catalog and that derived by Abdo et al. (2010a) is only 
about 20%, which is negligible for our study. We estimated the radii using the catalog of Galactic 
SNRs by Green (2014). 

The flux ratio between 0.3-10 keV and 0.1-100 GeV, Fx/Fq in Table 2, is in the range of 
~0.1-30. If 3FGL J1214.0—6236 is the gamma-ray counterpart of G298.6—0.0, our target has the 
faintest Fx/Fq ratio. Figure 3(a) shows the radius vs. Fx/Fq of MM SNRs with GeV gamma-rays. 
The flux ratio shows decline as the radii of SNRs is larger, or SNRs evolves older. The flux ratio of 
our target is well consistent with this trend. 

It is unclear why the flux ratio becomes smaller when MM SNRs evolve. One possibility 
is that the X-ray luminosity decreases faster due to the expansion and cooling (Gehrels & Williams 
1993), whereas the GeV gamma-ray luminosity does not decrease. In order to check this hypothesis, 
we have made plots of radius vs. 0.3-10 keV luminosity (Figure 3(b)) or 0.1-100 GeV luminosity 
(Figure 3(c)). This is consistent with our scenario; The X-ray luminosity become smaller when SNRs 
evolve (or radii become larger), whereas the gamma-ray luminosity is almost constant in the range 
of lO^^-lO^® erg s“^ This is also consistent with previous results in X-rays (Long 1983) and GeV 
gamma-rays (Brandt et al. 2015). The constant gamma-ray flux support our hypothesis that timescale 
for low-energy particles to escape these SNRs is long in comparison to the SNR evolution. This 
trend also implies that a part of Galactic GeV unID sources can be very evolved MM SNRs without 
significant detection in X-ray and radio bands. Note that radio surface brightness becomes smaller 
when the SNR evolves (Case & Bhattacharya 1998), which makes difficult to detect them together 
with large angular size. Further study of a large sample of GeV gamma-ray SNRs is needed to 
understand this tendancy. 
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Fig. 1. (a): The 843 MHz (contours) and 0.5-5 keV (gray scale) images of the G298.6-0.0 region. The X-ray image is smoothed with a Gaussian function 
with the kernei of 0.4 arcmin, Non-X-ray background (NXB) was not subtracted, and no vignetting correction was performed. The color scale is in logarifhmic 
scale. Red circle shows fhe error region of the Fermi source 3FGL Jf 214.0—6236. The coordinate is in J2000. The point source in the north-west is 
1RXS J121248.7-623027 (Fuhrmeister & Schmitt 2003). (b): Same image to (a) with source (bold) and background (thin) regions for the spectral analysis. 
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Fig. 2. Background-subtracted spectra of the source region. Black and red dosses represent FI and Bl spectra, respectively. The solid lines represent the 
best-fit VAPEc model. 
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Fig. 3. (a): Radius vs. flux ratio between 0.3-10 keV (absorption-corrected, Fx) and 0.1-100 GeV (Fa) for MM SNRs with GeV gamma-ray emission. 
Circles and the box represent the samples from literatures and from this work, respectively, in all panels, (b): Radius vs. 0.3-10 keV luminosity, (c): Radius 
vs. 0.1-100 GeV luminosity. 
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Table 2. X-ray and gamma-ray fluxes of MM SNRs with associations to GeV gamma-rays 


Sample Radius** Ex'" Fa‘^ Fx/Fq References 

W49B. 8 5.1 ±0.7 56" 1.8 31 (1) (2) (3) 

G349.7-I-0.2. . 11 3.0 ±0.4 4.1" 0.49 8.4 (4) (5) (3) 

3C 391 . 7.2 7.5 ±1.5 4.5 0.49 9.3 (6) (2) (3) 

W28. 1.8 13.1 ±2.6 1.1^ 3.6 0.30 (7) (8) (3) 

W44. 2.8 13.5 ±1.5 1.5^ 5.4 0.28 (1) (2) (3) 

G298.6-I-0.0. . 10 15 ±2 0.059 0.44 0.14 this work, (3) 


Note — (1) Claussen et al. (1997), (2) Kawasaki et al. (2005), (3) The Fermi- 
LAT Collaboration (2015), (4) Tian & Leahy (2014),(5) Yasumi et al. (2014), (6) 
Radhakrishnan et al. (1972), (7) Clemens (1985), (8) Rho & Borkowski (2002) 
“:Distance to the source in the unit of kpc. 

**: In units of pc. Errors indicate the radii along the major and minor axes, 
absorption-corrected 0.3-10 keV flux in units of 10“*^° erg cm“^s“^. 

0.1-100 GeV flux in units of 10“^° erg cm“^s“*^. 

": Only ISM component was considered. 

^: The X-ray flux was not extracted from the entire SNR but from the brightest 
central regions. 
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